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To test the strong-gravity regime in Randall- Sundrum braneworlds, we consider black holes bound 
to a brane. In a previous paper, we studied numerical solutions of localized black holes whose horizon 
radii are smaller than the AdS curvature radius. In this paper, we improve the numerical method 
and discuss properties of the six dimensional (6D) localized black holes whose horizon radii are larger 
than the AdS curvature radius. At a horizon temperature T « l/2n£, the thermodynamics of the 
localized black hole undergo a transition with its character changing from a 6D Schwarzschild black 
hole type to a 6D black string type. The specific heat of the localized black holes is negative, and the 
entropy is greater than or nearly equal to that of the 6D black strings with the same thermodynamic 
mass. The large localized black holes show flattened horizon geometries, and the intrinsic curvature 
of the horizon four-geometry becomes negative near the brane. Our results indicate that the recovery 
mechanism of lower-dimensional Einstein gravity on the brane works even in the presence of the 
black holes. 

PACS numbers: 04.50.+h, 04.70.Dy, 04.70.Bw, 11.25. Wx 



I. INTRODUCTION 



Much progress has been made in understanding the dynamics of gravity, brane and localized matter fields in the 
context of braneworlds. The braneworld models proposed by Randall and Sundrum (RS) Q,Q have an extra dimension 
^ i that is compactified via anti-de Sitter (AdS) space and branes, assuming that we live on a brane (or a domain wall) 
ON inside the AdS space. This compactification is called warped compactification. The RS models do not require a large 
extra dimension due to the warped compactification, contrary to other interesting braneworld models with large extra 
dimensions [1 

In the warped compactification, analyses of linear and second-order perturbations show that gravity induced on the 
brane by localized matter fields is consistent with usual four-dimensional Einstein gravity 0, IE Q>H,H> El fill IT2I 
IT3I Il4| . although the problem is essentially a five-dimensional one. For the strong-gravity reg ime in the RS models, 
the low energy (or derivative) expansion method would be a useful approach in some cases | l5lll6lll7| . A recent fully 
nonlinear numerical treatment of time-dependent braneworld models is also an interesting development to understand 
' , nonlinear brane dynamics |l8j . 

A pioneering approach to st rong ly gravitating systems in the RS model was developed by Wiseman in the study of 
relativistic stars on the brane |l9| . In this study, the full five-dimensional Einstein equations are solved numerically, 
r* ■ employing a conformal gauge method that reduces the Einstein equations to a suitable boundary value problem under 
. £h \ an assumption of static axisymmetry. A consequence of nonlinear effects in standard four-dimensional Einstein gravity 
is the existence of an upper mass limit at a fixed radius of relativistic stars. The existence of the upper mass limits is 
confirmed for small stars whose proper radii 1Z are less than the AdS curvature radius £. While the upper mass limits 
for large stars whose radii are 1Z < 3£ are not confirmed, all observed physical quantities indicate that an effective 
four-dimensional description by means of four-dimensional Einstein gravity can be applied far into the nonlinear 
regime. Then the so-called recovery mechanism of Einstein gravity on the brane would persist in the strong-gravity 
regime. 

The basic idea of the conformal gauge method is a key to exploring the nonlinear regime of gravity with extra 
dimensions. Indeed, it has been applied for nonuniform black strings [20L l2l| and Kaluza-Klein black holes p^. l23|. 
As is emphasized in Ref. p(j| . the problem of nonuniform black strings is a clean and simple application of the 
conformal gauge method compared with the other applications, due to a technical issue of numerics associated with 
the coordinate system at the symmetry axis (see, e.g., Refs. |24l l25| for a similar problem in numerical relativity, 
although the technical issue is intrinsic to the conformal gauge). 

The method is also applied to a study of black holes in the RS braneworld model [2|1 [2?J • Physically realistic black 
holes are thought to localize on the brane, as matter fields localize there [2^. Static axisymmetric localized black 
holes, for example, would be formed as a result of gravitational collapse on the brane, settling down to a static state 
at late times. There are many discussions about the localized black holes psl Ejj I30L l3lL 133. 133 . 134L I3EL l3tj l37l l38l I3H 
lifl liH H3. 143) (recent discussions are 0, ^5|, E3 ) ■ However, physically acceptable black hole solutions have not 
been found so far, except for exact solutions in the four-dimensional RS braneworld found by Emparan, Horowitz, and 
Myers (EHM) 0> HH ■ An interes ting arg ument motivated by a no-go theorem and the AdS / CFT correspondence 
in the RS braneworld (e.g. Refs. |5lll52j ) is that large localized static black holes may not exist |53l l54j . 
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In our previous papers, we considered a localized black holes in the five-dimensional RS single brane model, and 
presented some examples of numerical solutions H^,^. However, we could not succeed in constructing localized 
black holes that are larger than the AdS radius, because the numerical scheme gave no convergent solution for them. 
To overcome the difficulty and to obtain an understanding of localized black holes, we need further development of 
the method. A key method had been developed through the study of Kaluza-Klein black holes (see also Ref. 

In this paper, we apply the improved numerics and discuss numerical solutions of static localized black holes. 
We perform the analysis in the six-dimensional RS model. This is a technical reason that numerical calculations in 
higher dimensions have an advantage: The metric in higher dimensions decay faster than that in the original five- 
dimensional braneworld model, and which makes numerical calculations more tractable |20j| . We think that there is no 
qualitative difference between localized black holes in five-dimensions and those in six-dimensions. Indeed, as we will 
see, thermodynamic properties of the small localized black holes in five-dimensions found in Refs. |2a. |27| persist in 
six-dimensions. While the EHM solutions are known in four-dimensions, the difference between four-dimensions and 
dimensions more than four is very large. For example, the EHM black holes have a maximum mass, corresponding to 
the asymptotic deficit angle on the brane. 

We are interesting in qualitative properties of the localized black holes, for example, the thermodynamic property, 
horizon geometry and gravity induced on the brane. In particular, it is an interesting question whether the lower- 
dimensional Einstein gravity is recovered on the brane even in the presence of a localized black hole. We construct 
the static localized black holes by means of the numerical method and discuss the fundamental properties of the black 
holes. 

In Section[n] we begin by briefly reviewing the conformal gauge method, and then in Section [IIII we present numer- 
ical solutions, showing the performance of the method. Properties of the localized black holes and their implications 
are discussed in Section liVl Section[V]is devoted to a summary. Due to the technical nature of Section ITTT1 the reader 
may prefer to skip straight to Section llVI 



II. SIX-DIMENSIONAL BRANEWORLD BLACK HOLES 



A. Conformal gauge method 

Following previous work [2(||22j], we adopt the conformal gauge method to construct localized black holes. We begin 
with a general ansatz for static axisymmetric black hole solutions to Einstein's equations with a bulk cosmological 
constant in D dimensions: 

ds 2 = ^ (-e 2a dt 2 + e 2 (\dr 2 + dz 2 ) + r 2 e 2 W D _ 3 ) , (I) 

where the metric functions are arbitrary functions of r and z. The bulk cosmological constant is A — 
— (D — 1)(D — 2)/2£ 2 , and tension of the branes is a = 2(D — 2)/kd£. Here kd — 8ttGd is £>-dimensional New- 
ton's constant and I is the AdS curvature radius. Israel's junction condition at the brane in /^-dimensions gives 
K fJl v\ z= i+ = —j^v/i , where K^ v is the extrinsic curvature of the brane at z — t with the induced metric j^. 

Now we consider a six-dimensional (6D) single brane model. The metric form has residual gauge degrees of freedom 
of conformal transformations in the two-dimensional space {r,z}, and then we can use these conformal degrees of 
freedom to transform the location of the event horizon to be at a constant radius 



Here we have introduced polar coordinates, 



P = Po- (2) 



r = psmx, 

z = 1 + p cos x ■ (3) 



Note that the positive tension brane can be placed at z = £. We introduce 

£ = X 2 . (4) 

This variable is used in Ref. j2(J to gain numerical stability near the axis. However, even with this modification, the 
algorithm is still unstable and encounters a lack of convergence when we increase the radius of black hole. As we 
discuss later, this numerical problem is resolved in Ref. |22j |. altering a method originally used in Ref. 0]. An a i m 
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in using the variable Q is to increase a numerical resolution near the brane. All numerical calculations presented in 
this paper are performed by using the variable (0J . 

If the localized black hole is sufficiently small compared to the curvature radius in the bulk, the black hole would 
be approximated by the 6D Schwarzschild black hole because the tension of the brane and the bulk cosmological 
constant would not affect locally the black hole. This is in fact observed in the study of 5D braneworld black holes 
|27j . Thus we subtract the 6D Schwarzschild metric from the assumed metric. Rewriting the metric functions as |2^| 

a = A , 

(3 = B — C , 

7 = B + ^C , (5) 

and we decompose these into 

X (p,0=X s (p)+X(p,O. (6) 
where X = A,B, or C. Each X$ is defined by 

Cs(p) = 0. (7) 

When we take A = B = C = and £ — > oo (£/z — > 1), the 6D Schwarzschild solution in isotropic coordinates is 
recovered. 

Black hole solutions are parametrized only by a parameter 2 

Po , Q S 

Po = — ■ (8) 

As we will see later (Fig. EJ), the proper circumferential radius of the horizon on the brane is approximately equal to 
the parameter po- Black holes with po < I are small compared with the AdS curvature radius, whereas black holes 
with po > I are large. 



B S (P) = -log(l + 9) , 



B. Einstein equations 



The Einstein equations yield five nontrivial equations. Three of them yield elliptic equations for the respective 
metric components X = {A.B.C}, 

(d 2 r +d 2 z )X = S X , (9) 

where the source Sx depends nonlinear ly on X and its derivatives. The remaining two equations give constraint 
equations, which satisfy Cauchy-Riemann relations. The Cauchy-Riemann relations are obtained from the nontrivial 
components of the Bianchi identities VG^^ = 0, assuming that the equations for A, B, and C are satisfied. The 
Cauchy-Riemann relations in terms of are 

A/ ^<%$ + dp* = 0, (10) 



P 



a ? *-a p $ = o, (ii) 



1 The introduction of the variable £ is not essential in our numerics, but the increased resolution makes numerics work better, compared 
with calculations using the variable \- Indeed, this allows us to find large r bla ck holes. 

2 In our numerical calculations, we take po = 1 and change I to specify po l26l . Because the variation of I keeping po fixed corresponds 
to rescaling the length scale, we must be careful to compare physical quantities of a black hole to those of a black hole with different 
po. Thus in this paper we plot physical quantities in dimensionless form. We also note that physical size of the grid resolution (d^, dp) 
changes for each calculation due to the rescaling of the length scale, even if we fix the grid resolution. 
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where $ := 2^-AetgG^ and := py/lyj -&e\.g(G p p - G\). * and $ are called weighted constraint equations. Since 
\1/ and <& satisfy Cauchy-Riemann relations, each of them satisfies the Laplace equation. Hence if ^ is satisfied on all 
boundaries and provided that $ vanishes at any one point, the two weighted constraints must vanish in all places as 
long as the elliptic equations for A, B, and C are satisfied. 



C. Boundary conditions 

To solve the elliptic equations © numerically, we need boundary conditions for them. The boundary conditions 
are determined by physical requirements. Let us first consider the boundary conditions on the event horizon. The 
horizon is a regular singular point of the elliptic equations for fixed X- The regularity at the horizon demands that 
singular solutions should be absent. Thus a power expansion in p of the elliptic equations and the constraint equation 
{G p p — G|) = results in boundary conditions at the horizon for regularity: 

d p A = d p B = ^, d p C = 0. (p = p ) (12) 

z 

In addition to these conditions, the regularity of the constraint equation G^ = gives a condition 

d x (A-B + C)=0. (p = po) (13) 

This additional requirement is physically a condition that temperature of a black hole is constant on the horizon (see 
Eq. JUJ). We do not use this boundary condition for the elliptic equations, but we use it to integrate the constraint 
equation <E> (see Sec. IIII A|) . 

Boundary conditions on the brane are derived from the Israel's junction condition. They are given by 



d x A = d x B = -^ 



(x = tt/2). 

d x C = 0. (14) 



On the axis r = 0, the absence of conical singularity requires 

C(p,0)=0. (15) 

This condition is also required from the regularity of the equation for C on the axis. Furthermore, the axial symmetry 
requires that the r derivative of the metric functions vanishes on the axis. Thus we have 

d x A, d x B, d x C = 0. ( X = 0) (16) 

Let us consider the asymptotic boundary conditions. The background spacetime of the RS braneworld is the AdS 
spacetimc. It is well known that localized matter fields reproduce linearized lower-dimensional Einstein gravity on 
the brane in the weak gravitational regime, and the linear perturbations fall off asymptotically keeping the spacetime 
to be asymptotically AdS. Thus in the current problem, we could assume that the spacetime is asymptotically AdS. 
Therefore, the asymptotic boundary conditions are 

A,B,C = (p^oo) (17) 

at spatial infinity, although in actual numerics these conditions are imposed at a finite coordinate location. In Appendix 
IB1 we discuss asymptotic behavior of gravitational field induced by a point mass on the brane. The gravitational 
field in the asymptotic region is anisotropic in the polar coordinates [a, Q, |55| . While the asymptotic behavior given 
by the point mass could provide asymptotic boundary conditions for the localized black holes as Neumann boundary 
conditions, we have not been able to find a stable scheme involving the boundary conditions. However, already with 
our simple implementation it is possible to derive interesting results as we will see (see Sec. 1111 HI for tests of our 
numerical method) 



III. NUMERICAL SOLUTIONS 
A. Stabilization of numerical scheme 



In this section, we discuss the stability of our relaxation scheme of numerics. A problem of solving the elliptic 
equations @ by a relaxation scheme is an instability coming from the axis (r = 0). A term in the equation for C 
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FIG. 1: Metric functions and C — C2 for po = 0.5. 





<l*l> 
po = 0.5 1.25 


(1*1) 

0.5 1.25 


{\C-C 2 \) 

0.5 1.25 


32 x 128 


4.1 xlO^ 3.4 xlO -1 


1.9 xlO - * 4.0 xlO - * 


3.1 xlO~ a 7.1 xlO~ a 


64 x 256 


1.1 xlQ'' 2 7.8 xl0~ 2 


4.9 xlQ- A 1.6 xl(T a 


8.7 xlO~ 4 2.0 xlO~ a 


128 x 512 


3.5 xl0~ a 1.4 xl0~ 2 


1.6 xHT 3 3.7 xHT 3 


2.5 xlO~ 4 5.5 xlO~ 4 



TABLE I: Averaged violations of the absolute values of the weighted constraints and C — C2 for the three resolutions. The 
three resolutions converge consistently with a second-order scaling. The averaged violations are calculated for a fixed domain 
1 < p/p < 5. 



that has a factor 1/r 2 gives severe contribution near the axis, and it invokes numerical instability |2fij . Actually if we 
turn off the term during a calculation, then the calculation becomes very stable and converges to a solution, although 
the solution is not a solution of Einstein's equations. This issue of instability prevents us from finding large localized 
black holes in our previous work. 

We deal with this term as in Ref. ■ A second derivative term for C in the constraint equation <I> oc is simply 
C tP £, and thus we can integrate it over the domain. A point is that this equation has no singular term as long as 
Eq. I|lfi(l is satisfied. Thus the solution for C integrated away from the axis has very good behavior near the axis. 
We call the solution C2 to distinguish it from C. Given the characteristics we need two initial data surfaces, one at 
constant £, and the other at constant p. Suitable surfaces are £ = and p = p$. For £ = we can set C2 — to 
satisfy the condition l|15|l . and this ensures the good quality of the C2 behavior near the axis. For the other initial 
surface p = p Q , we fix C2 using the condition (| 1 3|) that the horizon temperature is a constant. We then integrate C2 
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FIG. 2: Metric functions and C - C 2 for p = 1.25. 



1*1 



1*1 




r r 

FIG. 3: Absolute values of the weighted constraint equation \& for po = 0.5 (left) and po = 1.25 (right). 



from these two boundaries: 



C(po,Q = (B - A) - (B - A)\s =0 . 



(18) 



where F is the remaining terms in the constraint equation. We calculate the 1/r 2 term in the elliptic equation for C 
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5(1*1) 


6{\*\) 


ST 

n ^ 1 w 


SS 6 

n ^ 1 w 

U.U, 1 


U.O, J- . £u 


Pmax/PO = 6.0 


0.673 0.369 


12.9 1.32 


-0.014 -0.029 


0.082 0.133 


0.071 0.126 


7.2 


-0.020 0.319 


3.12 0.594 


-0.009 -0.017 


0.049 0.081 


0.043 0.077 


9.7 


0.016 0.145 


-0.104 0.164 


-0.003 -0.006 


0.016 0.027 


0.014 0.026 


11.0 


0.011 0.060 


0.000 0.155 


-0.001 -0.003 


0.007 0.012 


0.006 0.011 



TABLE II: Sensitivity of the solutions to the finite position of the asymptotic boundary p ma .x- We compare the temperature, 
entropy, and weighted constraints for solutions with po — 0.5 and 1.25. The variation of the quantities is defined as ST = 
^"/^(pmax =12.2) — 1, dividing by T for p ma x = 12.2po- The averaged absolute values of the weighted constraints, (l^l) and 
(l^l), are calculated for a fixed domain, 1 < p/ po < 5. The constraints appear to be even better satisfied as pmax increases. 
Moreover, thermodynamic quantities hardly change, indicating our choice of p ma x is sufficiently large. All solutions are at the 
same resolution as for 64 x 256 with p ma x = Upo- 



using this function G2, rather than G. Then the term in question is replaced by 

^(l-^)-^(l-e^). (19) 

This prescription dramatically improves stability of numerics, and we are free of the instability coming from the 
axis. Furthermore, we use condition l|13|) in our scheme, the variation of the horizon temperature is vanishingly small 
(< 10~ 2 %), and the zeroth law of black hole thermodynamics is well satisfied. 

B. Performance of the method 

The formulations and numerical scheme described in the previous sections allow us to solve Einstein's equations 
as a boundary value problem. To solve the problem numerically, we use a relaxation method with a second-order 
discretization scheme. In this section, we present some results of numerical calculations and discuss qualitative aspect 
of numerical solutions. 

An illustration of metric functions for po — 0.5 and po = 1.25 are shown in Figs. and [21 The former case is a 
typical for small localized black holes (po < 1) and the latter case is for large localized black holes (po > 1). As we 
see, the global behavior of metric functions hardly change for the two cases, but the peak values of metric functions 
increase, indicating that nonlinear effects become much more important for large black holes. It is remarkable that 
C — C2 is very small, and then the integration method works consistently. We have performed the calculations for 
three different resolutions, 32 x 128, 64 x 256, and 128 x 512, setting d(, — dp (see Tablc[|J. The finite position of the 
asymptotic boundary is at 

Pmax = ll/Oo- (20) 

In our simple implementation, largest horizon radius of a black hole that we could find depends on the resolution of 
the calculation, and then it is limited by computation time. The largest black hole that we found within a moderate 
computational time is po — 2 with the grid resolution 128 x 512. The localized black holes are systematically 
constructed for the parameter region po = 0.002 ~ 2. 

Figure shows absolute errors of the weighted constraint equation "J cx (G p p — G|). Because $ ex is essentially 
the same as \I/ due to the Cauchy-Riemann relations and we have confirmed the consistency of G^ by showing C — C2 

in Figs. ^] and |21 we suppress figures of $. While the corresponding Einstein equation G p p — G| = is well satisfied 
yielding vanishingly small values, peaks of the errors of the weighted constraints appear in the asymptotic region 
near the brane. The peaks are originated from the huge amplification of the tiny violations of Einstein's equations 
due to the weight function p^—g that is used to define \F The weight function is proportional to ~ p 5 near the 
brane. Because the tiny violation of Einstein's equations comes from the asymptotic boundary, we will need to avoid 
using the naive asymptotic boundary conditions 1|17|) . As we have discussed, the naive Dirichlet conditions could be 
replaced by appropriate Neumann conditions. The prescription will reduce such errors from finite boundary effects, if 
we can implement it without any numerical instability. Nevertheless, as p m ax increases, the approximation using the 
Dirichlet boundary conditions (|17H becomes much better, and then violations of the weighted constraints decrease. 
This is observed in our numerics (Tablc^J. The averaged absolute errors of the weighted constraints in a fixed domain, 
(|$|) and (l^l), decrease as p max increases. Thus, we do not care about the violations of the weighted constraints in 
the asymptotic region in the present numerics. In Table ITT1 we also show that errors from finite boundary effects in 
thermodynamic quantities are < 1% for Eq. H20fl . and then the location of the boundary is sufficiently far. 
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IV. PROPERTIES OF THE LOCALIZED BLACK HOLES 



A. Thermodynamic quantities 

We have demonstrated the numerical method and shown numerical examples of the localized black holes. The 
largest localized black hole that we found is po = 2. In this section, we try to understand the properties of the 
localized black holes. Let us first consider thermodynamic quantities of the localized black holes. Thermodynamic 
quantities of interest are the horizon temperature T and the entropy Sq of the black holes, which are given by 



S 6 = 



1 

1 3 



2pt J dn 3 J* /2 d X (^V°+ c °sin 3 X 



(21) 



In Fig. 01 we plot a thermodynamic relation between these quantities normalized by the AdS radius I. The factor 
Gq/£ 4 is introduced to obtain correctly rescaled quantity (see footnote in Sec. Ill A|) . Three resolutions are used to 
generate the figures. One sees that they are consistent and compatible with second-order scaling. In the figure, we 
have also plotted the same thermodynamic relations for the 6D Schwarzschild black hole and the 6D black string, for 
comparison (see Appendix 0. One finds that for sufficiently small horizon radius (T ^> t ), the thermodynamic 
relation of the localized black hole is a 6D Schwarzschild black hole type. This behavior is confirmed for the 5D 
localized black holes j23. On the other hand, as the horizon radius increases (T <C i^ 1 ), the thermodynamic relation 
becomes a 6D black string type. The transition point at which the localized black hole changes its thermodynamic 
character is given by a temperature 

r-jL. (23) 

As could be expected, this temperature is given by the temperature of the localized black hole with po — 1, which is 
T\p 0= i 0.9/2ir£. It is interesting that the horizon temperature l|23ll is the same as that of a "critical" black string 
whose horizon radius is = I. Furthermore, the temperature reminds us of the Hawking-Page phase transition 
of the Schwarzschild-anti-de Sitter black hole [56j. 

We consider the specific heat, TdSg/dT, of the localized black hole. It is immediately clear from Fig. Q] that 
the specific heat is negative. Because the temperature T and the entropy 5g are well approximated by the 6D 
Schwarzschild black hole and the 6D black string in the respective regions (T < I and T > i) , the specific heat of the 
localized black hole is also well fitted by them. 

There is an interesting quantity that characterizes the localized black hole by a geometrical quantity on the brane 
|27|. It is black hole "entropy" measured by an observer living on the brane. We define the five-dimensional entropy 
by using proper circumference at z = £ as 

S 5 = ^[2n 2 ple 3B ^ c °] z=e . (24) 

In Fig. we display a relation between S§ and T taking an appropriate non-dimensional combination. In the figure, 
the same thermodynamic relations for the 6D black string and the 6D Schwarzschild black hole are also plotted for 
comparison. The figure shows that the five-dimensional entropy S§ deviates from that of the 6D Schwarzschild black 
holes as the temperature decreases. Then the entropy 6>5 tends to behave like that of the 6D black string. Because 
of limited resolution and computing power, we could not confirm whether S§ converges to the value of the 6D black 
string in the limit of large horizon radius (T — > 0). However, an extrapolation of our results for S§ and Sq suggests 
that the thermodynamic relations of the localized black hole converge to those of the 6D black string in the limit. 
Recall that the gravity induced on the brane by the 6D black string is the 5D Schwarzschild black hole. Therefore, 
if the convergence is true, it means that the gravity induced on the brane by the large localized black holes is well 
described by the 5D Einstein gravity. Thus, the above results imply that the recovery mechanism works even in the 
presence of the black hole, which is not confined on the brane but localizes there, extending into the bulk. 

We give a comment on a possible deviation from the thermodynamic character of the black string. For the weak 
gravity regime in the original 5D braneworld model, the relative order of the Kaluza-Klein corrections to Newton's 
law is 0(£ 2 /r 2 ) where r is a coordinate radius. Similar Kaluza-Klein corrections will appear in the gravity 

induced on the brane by a black hole. However, the correction will be very small for massive black holes, and then 
numerical calculations can probe only the leading behavior. Therefore, even if numerical calculations show a good 
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FIG. 4: Thermodynamic quantities of the localized black holes. The left figure shows a thermodynamic relation for the 
temperature T and the entropy <5>6, which are normalized by the AdS curvature radius i. The right figure shows an appropriate 
combination of the 5D entropy S5 and the temperature T as a function of the temperature. The thermodynamic relations 
for the 6D Schwarzschild black hole (dot-dashed line) and the 6D black string (solid line) are also plotted for a reference. To 
generate the figures, three resolutions are used, 32 x 128 (cross), 64 x 256 (circle), and 128 x 512 (plus). 




-8 -6 -4 -2 
\n[MG 6 A 3 ] 



FIG. 5: The relation between the thermodynamic mass and the entropy for the localized black holes (solid curve). The dashed 
and the dot-dashed lines represent the relations for the 6D Schwarzschild black hole and the 6D black string, respectively. 



agreement between the thermodynamic character of the localized black hole and that of the black string, it does not 
mean that no correction to the 5D Schwarzschild black hole appears on the brane. 

Mass of a black hole is an important quantity to characterize the black hole. If one assumes the first law of 
thermodynamics psl l4fll | , one can calculate a thermodynamic mass of a localized black hole. Because the localized 
black hole has the thermodynamic character analogous to those of the 6D Schwarzschild black hole and the 6D 
black string, the thermodynamic mass obtained by integrating the first law also agrees with them in the respective 
regime (Fig[5J). The mass of large black holes is then approximated by that of the 6D black strings, neglecting small 
corrections coming from the integration over small radii. It is given by 

1 / Q \ 2 / 3 / / \ V3 

The black hole entropy Sq is greater than or nearly equal to that of the black string with the same mass. This is 
contrary to the argument for the (BTZ-like) localized black holes and the BTZ black strings in Ref. • 
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FIG. 6: Ratio of the equatorial to polar circumferential radii for the localized black holes (po < 2). Three resolutions are used 
[32 x 128 (cross), 64 x 256 (circle), and 128 x 512 (plus)]. The solid curve represents the equatorial radius R eq with respect to 
Po- 




FIG. 7: Intrinsic 4D curvatures of the horizon four-geometry for pa = 0.5, 1, 1.4. The curve showing the largest positive 
curvature on the axis £ = corresponds to po = 0.5. As the radius increases, the curvature on the axis decreases. While the 
curvature remains positive for po < 0.5, the curvature becomes negative near the brane for po > 0.5. 



B. Horizon geometry and curvatures 

A localized black hole is expected to have a flattened horizon geometry, and it will tend to flatten much more as 
its horizon radius increases |26i l27l Esj . Let us define polar and equatorial circumferential radii of the horizon by 



R eq = re Ba+ i Co \ z= 



4p 



t/2 



o B — Co 



(26) 



^polar — t . 

Z7T 

Figure shows R e q/Rpoiar and R eq . The small black holes are almost spherical, R eq m R po iar, but the large black 
holes show flattened geometry, R eq > R po iar- 

The intrinsic curvature of the horizon geometry gives us another information. The horizon four-geometry is 



dS 2 = — 



„2(B -C ) p 2 dx 2 + r 2 e 2B 0+ AC Q /3 dQ 2 



(27) 



The 4D intrinsic curvature R^' of the four-geometry is shown in Fig. 0for several black holes. The curvature becomes 



11 




FIG. 8: The four-geometry of the event horizon embedded in the Euclidean space HBj; based on our numerical solutions. Two 
rotational degrees of freedom are suppressed. The size of the localized black holes in the left and right figures are po = 0.1 and 
0.5, respectively. The corresponding equatorial circumferential radii are R eq /£ = 0.15 and 0.63, respectively. The solid curves 
show the location of the orbifold brane where the Zi. symmetry is imposed in the original space. 



negative near the brane for the black holes with po > 0.5. The appearance of the negative spatial curvature near the 
brane is also confirmed for the EHM solution and the 6D black string (Appendix |A"|) . Because of this negative spatial 
curvature, embeddings of the spatial horizon geometry into a five-dimensional Euclidean space 

dE 2 = dX 2 + dY 2 + Y 2 dn 2 3 (28) 

are not possible for all black holes, but only possible for small black holes (see, e.g., Ref. The cylindrical and 

radial coordinates of the horizon are given by 



X bh ( X ) = 




Y bh ( X ) = -e B «+i c «. (29) 

z 

Figure[5]shows examples of the horizon geometry embedded in the Euclidean space. Because the cylindrical coordinate 
of the horizon on the axis does not coincide with the polar circumferential radius R po iar, the embeddings in the 
Euclidean space do not conserve the ratio R eq / Rpoiar- To illustrate the horizon geometry, we rescale the cylindrical 
coordinate X for each black hole to conserve the ratio. The rescaled cylindrical coordinate is denoted by Xi (Fig. [SJ|. 
This illustration allows us to have an intuitive look of the horizon geometries for different horizon radii. 

While we have compared the localized black hole to the black string, the black string has a curvature singularity 
at the AdS horizon as well as the ordinary singularity at the center of the event horizon (Appendix [SJ 2S| . Thus 
the solution has been considered as an unphysical solution. Do the localized black holes have a naked singularity 
somewhere? To test the regularity of the spacetime, we calculate the square of Riemann tensor that is one of curvature 
scalars, 

K (6) ._ R n»\p Rfu/Xp . (30) 

It is hard to numerically calculate this kind of curvature scalars with enough accuracy near the axis and the horizon. 
In particular near the axis (r = 0), we encounter severe cancellations between terms that have coefficients proportional 
to 1/V 4 , 1/r 3 , • • • . Nevertheless, we can confirm that the curvature K^ 6 ' is finite with good behavior. Then as far as 
the curvature is concerned, localized black holes have no naked singularity, at least in the computational domain. 
This supports that our localized black holes are physically acceptable black hole solutions. In Fig. EH we illustrate 
the curvature on the brane. It goes rapidly to a constant value given by the AdSg, 



K^ = ™, (AdS 6 ) 



(31) 
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FIG. 9: 6D and 5D squares of Riemann curvature, \n(£ 4 K {6) - 60) (asterisk) and ln{£ 4 K^) (solid curve), on the brane for 
po = 0.5, 1.25. The dot-dashed curves show fits to the asymptotic behavior 1331 1. 



and then it begins to fluctuate due to roundoff error. Thus in Fig. we show for an appropriate range of proper 
circumferential radius R eq on the brane. We have furthermore calculated the square of the 5D Riemann tensor 

by using the induced metric on the brane. This curvature is a good measure for an observer on the brane. The scalar 
curvature for the 5D Schwarzschild black hole is 

K {5) K^r, (33) 

at a distance from the horizon [see Eq. l|A3jl ] . On the other hand, the curvature K 1 - 5 ' for the 6D Schwarzschild black 
hole is proportional to ~ We could then compare these behaviors with that of the curvature for the localized 

black holes. We plot K ^ of the localized black holes with po = 0.5 and 1.25 in Fig. [5] The dot-dashed curves are 
fits for the form (|33[1 in an appropriate asymptotic region. Although the localized black hole with po — 1.25 is not so 
large, we observe that the 5D curvature for po = 1.25 shows better agreement with Eq. 1(33(1 . compared to that for the 
smaller black hole. We note that the fit to the behavior 1)33(1 is better than a fit to the behavior oc However 
we refrain from going into a discussion of such asymptotic behavior, because the asymptotic boundary conditions 
adopted in this paper are too naive to discuss it. 



V. CONCLUSION 



No realistic black holes localized on the brane have been found so far in the warped compactification, except for the 
EHM solutions in the four-dimensional braneworlds. Black holes are important physical objects to understand and 
to test strong gravitational phenomena in the warped compactification. We have explored this problem by means of 
numerics, constructing black holes localized on the brane. We have employed the conformal gauge method that has 
been developed to study relativistic stars on the brane, nonuniform black strings and Kaluza-Klein black holes. In 
our previous work, we could not succeed in constructing large black holes, but could construct only the small black 
holes (po < 1/5). However, improved numerics allow us to construct the large localized black holes with po > 1, and 
then we have studied the properties of the localized black holes. 

The largest size of the black holes that could be relaxed in a reasonable computational time is po = 2 in this paper. 
The parameter po is approximately equal to the proper circumferential radius on the brane. The largest size depends 
on the resolution of the calculation, and thus it is limited by the computational time. The situation would be similar 
to the study of the large relativistic stars [lflj . 

The horizon geometry of the localized black holes is flattened in the bulk. The ratio of the equatorial to polar 
circumferential radii shows that the localized black hole becomes much more flattened as its horizon radius increases 
(Fig. |SJ) . We have pointed out that the intrinsic curvature of the horizon four-geometry becomes negative near the 
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brane for the black holes with po > 0.5. Then the global shape of the flattened horizon changes as its horizon radius 
increases. 

The thermodynamic properties of the localized black holes have been discussed by investigating the entropy, tem- 
perature, and specific heat of the black holes. For the small black holes (po < 1) the thermodynamic relation between 
the entropy and the temperature is the 6D Schwarzschild black hole type, whereas for the large black holes (po > 1) 
the thermodynamic relation is the 6D black string type (Fig. 0J. The transition occurs at T l/2w£, corresponding 
to po ~ 1- The 5D entropy measured by an observer on the brane also shows the transition: the physical quantities of 
the large localized black holes are approximately described by the 6D black strings. Thus the large localized black hole 
looks like the 6D black string with its ends capped off. This transition of thermodynamic relations might be analogous 
to an expected transition between a sequence of (static) black string solutions and that of black hole solutions in the 
Kaluza-Klein compactification (5^, although the transition is questioned by a recent study (Ref. f22\). 

Recall that the induced gravity of the 6D black string is given by the 5D Schwarzschild black hole. Then the above 
results indicate that the gravity induced on the brane by the large localized black hole is approximated by the 5D 
Einstein gravity. Such a recovery mechanism of lower-dimensional Einstein gravity on the brane is proved in the weak 
gravity regime by means of linear and second-order perturbations 

0,112, 

and it is partially confirmed for the strong 
gravity regime by the numerical study |19|. Our result also suggests that the recovery mechanism works even in the 
presence of the black hole that is bound to the brane. 

Let us discuss the relation between our results and the argument of the classical black hole evaporation conjecture 
[53L IH^ . The argument claims that large localized static black holes may not exist. On one hand, the physical 
quantities calculated in this paper have not indicated a critical point at which a sequence of solutions suddenly 
disappears or some kind of singular behavior appears. Thus it seems to shed doubts on the conjecture. On the other 
hand, however, the largest localized black hole in this paper is po = 2, and the black hole is not sufficiently large. It 
would be still small compared with the black holes supposed in the conjecture. It is certainly true that we cannot 
provide a conclusive disproof (or proof) of the conjecture by means of numerics, even if we construct a very large black 
hole with, e.g., po = 100. In Ref. [23], some possible consistent scenarios are discussed. However, we expect that the 
sequence of our localized black holes represents the unique stable vacuum solution without the naked singularity that 
describes a nonrotating uncharged black hole bound to the brane. Needless to say, it is preferable to have an analytic 
solution of a localized black hole, and thus much progress in the study of braneworld black holes is required. 

An interesting application of the conformal gauge method is a study of black holes in the RS two-branes model. 
Because the localized black holes in the model are similar to black holes in the Kaluza-Klein compactification, the 
method in Ref. can be directly applied. Effects of the warped compactification, the second brane, and radius 
stabilization mechanism on the black holes would result in different phenomena observed in this paper. Charged 
localized black holes are also interesting. Introducing gauge fields on the brane is not difficult in our approach. We 
hope to report soon our results toward these directions. 
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APPENDIX A: 6D BLACK STRING 

The black string solution in the single brane model has a singularity at the AdS horizon (z — > 00) 28] . However, 
the solution is still useful as a reference. The 6D black string is given by 

C 2 

ds 2 = — [-Vdt 2 + y- x dr 2 + r 2 dn 2 3 + dz 2 } , (Al) 



where V — 1 — r 2 /r . The brane is at z — t. The square of the Riemann tensor is 



60+^ 



(A2) 



which diverges at the AdS horizon (z — > 00), as well as at the black string singularity at r = 0. The 5D curvature 
scalar (l' ; >21 on the brane is given by 

if( 5 ) = 72r 4 h /r 8 . (A3) 
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The intrinsic curvature R^ of the horizon four-geometry is 

6 



R^ 



1(2 



(A4) 



This is not always positive, but becomes negative for z < \/2rh- 
The temperature and entropy of the black string are 



T = 



1 



1 

4G^ 



(27T) 2 



(A5) 



Assuming the first law, thermodynamic mass is given by 



2/ .3 



1/3 



APPENDIX B: ASYMPTOTIC BOUNDARY CONDITIONS 
1. Linear perturbations 

The asymptotic behavior of black hole characterizes the black hole itself, and then it is important. To impose 
physically reasonable asymptotic boundary conditions, we need to solve linear perturbations of the black hole geometry 
in the asymptotic region. Here we assume that the asymptotic behavior of black hole geometry is equivalent to that 
of gravitational field induced by a point mass on the brane. Then we consider solutions to the linear perturbations of 
weak gravitational field in the conformal gauge. Linearized equations in the conformal gauge do not yield decoupled 
equations, and then it is difficult to solve it. However, using the transverse-traceless gauge, one finds that a master 
equation exists and it is solved using a Green's function. The discussion in this appendix is based on an argument in 
Ref. 

Perturbations in the transverse-traceless (TT) gauge are given by 



dsf = — [-e z± 'dr + , 

where F and D are determined by the TT conditions 

1 



2F,,2 , „2H,2 . 2 2D 



dr 2 + r 2 e 2D dQ + dz 2 ] 



F =-Z3 d r(r 4 H), D=—d r (r 3 H). 



The linearized equation for H is given by 



z 2 



a 2 + -d r + d 2 - ~d z 

r ' z 



H . 



(Bl) 



(B2) 



(B3) 



Comparing these coordinates to those in the conformal gauge, coordinate transformations to bring the metric into 
the conformal gauge are r' TT = r + /, z' TT = z + g with g^ z = f^ r + H, g, r = —f tZ which yield Poisson equations for / 
and g 



Ag = H. , Af = -H. 



(B4) 



The linearized equation (|B3|I does not change by the coordinate transformations. Then the coordinate transformed 
metric components are given in terms of H, /, and g: 



A 
Bo 



fJ 



F , 



3/ 9 2^ 3^ 
■z- -- + 1-0*9+ -D. 
or z 5 5 



r 3 ^ 3 ^ 3 n 



(B5) 
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To solve the coordinate transformations (|B4ll for / and g, we need to specify the boundary conditions. Boundary 
conditions on the brane are derived from the brane matching condition that is expanded to linear order, 

d 2 p g + -d p g = 0, 

3d p g = H x , (B6) 

and we also have 39 x / = d x H . On the axis, we impose 

/ = 0. (B7) 

The asymptotic behavior l|B5(l of gravitational field induced by a point particle can be used as asymptotic boundary 
conditions of black hole geometry. The boundary conditions are given by Neumann conditions, eliminating an overall 
coefficient that corresponds to a mass of a black hole. To give the boundary conditions i|B5|) . we must first solve 
Eq. (|B3|) and then determine the coordinate transformations by solving the Poisson equations l|B4|) . In the next 
subsection, we discuss a solution to the linearized equation (|B3|I and the coordinate transformations. 



2. Solution to linear perturbations 

Linearized gravity in the RS background has been discussed several times in the literature. A treatment of it in 
general dimensions is discussed in Ref. |55j| . In this reference, metric perturbations in the transverse-traceless gauge 
are solved by means of a Green's function. While the solution is general, analytic estimations of H and the coordinate 
transformations become somewhat messy. Thus to evaluate the leading order contributions, we consider another 
method that is rather simple. 

A solution of (|B3|) is constructed as 

fo / z\ 5/2 / / \ 2 f°° r . n 

H = l ^{j) (r) Jo dk [( M ) 3/2h ( k )Mkr)K 5/2 (kz)\ , (B8) 

where J 2 and K§/ 2 are the Bessel function and modified Bessel function, respectively, and h(k) is an arbitrary function. 
Because we are interesting in asymptotic behavior (r>£orz>f), only long- wave modes give dominant contribution 
in the integration. Let us consider h{k) = const. = ho- Then we obtain an approximate solution, 

_ hof_ /3 z 2 z 5 \ 

~ r 2 I 2 _ r 2 r 2 (r 2 + z 2 ) 3 / 2 J ' 



Particular solutions of Eq. (|B4|I for the asymptotic behavior (|B9ll can be found easily by using polar coordinates. 
They are 

h £ 2 (-8r 4 - 4r 2 z 2 + z 4 - zy/r 2 + z 2 (3r 2 + z 2 )) 
3r 2 {r 2 + z 2 )i 

To obtain solutions that satisfy required boundary conditions, we need to add homogeneous solutions. Then general 
solutions are given by 



g(r, z) = g p + dqF(q) cos(qr / t)e- qz '\ 
Jo 

poo 

f(r, z) = f p - dqF(q) 8m(qr/t)e-<>' e . (Bll) 
Jo 

From the boundary condition 3d r g = —rd z H on the brane, we obtain 

= ^ + *h£* [2K{ ) _ K{ )h (B12) 

3 6TT 

On the brane (z =f,r><), the leading order contributions are evaluated as 

£ 3 h 

9 * 
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-0(log(r/*)). 



(B13) 



Near the axis (z ^> £, r <C £), we have 



5 
/ 



O 

rhp 

z 



\ogz 



O 



O 



\ogz 



(B14) 



The coordinate transformation i|Bll|) and the asymptotic behavior (|B9|I (or the solution in Ref. [55]) with Eq. I|B5 
can be used as the asymptotic boundary condition of localized black holes. 



[1] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999), hep-ph/9905221. 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999), hep-th/9906064. 

[3] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Lett. B429, 263 (1998), hep-ph/9803315. 

[4] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Lett. B436, 257 (1998), hep-ph/9804398. 

[5] T. Shiromizu, K.-i. Maeda, and M. Sasaki, Phys. Rev. D62, 024012 (2000), gr-qc/9910076. 

[6] M. Sasaki, T. Shiromizu, and K.-i. Maeda, Phys. Rev. D62, 024008 (2000), hep-th/9912233. 

[7] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 (2000), hep-th/9911055. 

[8] T. Tanaka and X. Montes, Nucl. Phys. B582, 259 (2000), hep-th/0001092. 

[9] S. Mukohyama, Phys. Rev. D65, 084036 (2002), hep-th/0112205. 
[10] S. Mukohyama and L. Kofman, Phys. Rev. D65, 124025 (2002), hep-th/0112115. 
[11] I. Giannakis and H.-c. Ren, Phys. Rev. D63, 024001 (2001), hep-th/0007053. 
[12] H. Kudoh and T. Tanaka, Phys. Rev. D64, 084022 (2001), hep-th/0104049. 
[13] H. Kudoh and T. Tanaka, Phys. Rev. D65, 104034 (2002), hep-th/0112013. 
[14] H. Kudoh and T. Tanaka, Phys. Rev. D67, 044011 (2003), hep-th/0205041. 
[15] T. Wiseman, Class. Quant. Grav. 19, 3083 (2002), hep-th/0201127. 
[16] S. Kanno and J. Soda, Phys. Rev. D66, 083506 (2002), hep-th/0207029. 
[17] T. Shiromizu and K. Koyama, Phy. Rev. D 67, 084022 (2003), hep-th/0210066. 
[18] J. Martin, G. N. Felder, A. V. Frolov, M. Peloso, and L. Kofman (2003), hep-th/0309001. 
[19] T. Wiseman, Phys. Rev. D65, 124007 (2002), hep-th/0111057. 
[20] T. Wiseman, Class. Quant. Grav. 20, 1137 (2003), hep-th/0209051. 
[21] T. Wiseman, Class. Quant. Grav. 20, 1177 (2003), hep-th/0211028. 
[22] H. Kudoh and T. Wiseman (2003), hep-th/0310104. 
[23] E. Sorkin, B. Kol, and T. Piran (2003), hep-th/0310096. 
[24] M. Alcubierre et al., Int. J. Mod. Phys. D10, 273 (2001), gr-qc/9908012. 
[25] M. Shibata, Prog. Theor. Phys. 104, 325 (2000), gr-qc/0007049. 

[26] H. Kudoh, T. Tanaka, and T. Nakamura, Phys. Rev. D68, 024035 (2003), gr-qc/0301089. 
[27] H. Kudoh, Prog. Theor. Phys. 110, 1059 (2004), hep-th/0306067. 

[28] A. Chamblin, S. W. Hawking, and H. S. Reall, Phys. Rev. D61, 065007 (2000), hep-th/9909205. 

[29] A. Chamblin, H. S. Reall, H.-a. Shinkai, and T. Shiromizu, Phys. Rev. D63, 064015 (2001), hep-th/0008177. 

[30] C. Cadeau and E. Woolgar, Class. Quant. Grav. 18, 527 (2001), gr-qc/0011029. 

[31] R. Casadio and L. Mazzacurati, Mod. Phys. Lett. A18, 651 (2003), gr-qc/0205129. 

[32] R. Casadio, A. Fabbri, and L. Mazzacurati, Phys. Rev. D65, 084040 (2002), gr-qc/0111072. 

[33] N. Dadhich, R. Maartens, P. Papadopoulos, and V. Rezania, Phys. Lett. B487, 1 (2000), hep-th/0003061. 

[34] J. Garriga and M. Sasaki, Phys. Rev. D62, 043523 (2000), hep-th/9912118. 

[35] I. Giannakis and H.-c. Ren, Phys. Rev. D63, 125017 (2001), hep-th/0010183. 

[36] M. S. Modgil, S. Panda, and G. Sengupta, Mod. Phys. Lett. A17, 1479 (2002), hep-th/0104122. 

[37] P. Kanti and K. Tamvakis, Phys. Rev. D65, 084010 (2002), hep-th/0110298. 

[38] G. Kofinas, E. Papantonopoulos, and V. Zamarias, Phys. Rev. D66, 104028 (2002), hep-th/0208207. 

[39] H. Kodama, Prog. Theor. Phys. 108, 253 (2002), gr-qc/0204042. 

[40] S. Shankaranarayanan and N. Dadhich (2003), gr-qc/0306111. 

[41] T. Tamaki, S. Kanno, and J. Soda (2003), hep-th/0307278. 

[42] S. I. Vacaru and E. Gaburov (2001), hep-th/0108065. 

[43] S. I. Vacaru (2001), hep-th/0110250. 

[44] C. Charmousis and R. Gregory, Class. Quant. Grav. 21, 527 (2004), gr-qc/0306069. 
[45] P.-J. De Smet, Class. Quant. Grav. 20, 2541 (2003), gr-qc/0302081. 
[46] K. T. P. Kanti, I. Olasagasti (2003), hep-th/0307201. 

[47] P. S. D. Karasik, C. Sahabandu and L. C. R. Wijewardhana (2003), gr-qc/0309076. 



17 



[48] R. Emparan, G. T. Horowitz, and R. C. Myers, JHEP 01, 007 (2000), hep-th/9911043. 

[49] R. Emparan, G. T. Horowitz, and R. C. Myers, JHEP 01, 021 (2000), hep-th/9912135. 

[50] M. Bruni, C. Germani, and R. Maartens, Phys. Rev. Lett. 87, 231302 (2001), gr-qc/0108013. 

[51] M. J. Duff and J. T. Liu, Phys. Rev. Lett. 85, 2052 (2000), hep-th/0003237. 

[52] S. W. Hawking, T. Hertog, and H. S. Reall, Phys. Rev. D62, 043501 (2000), hep-th/0003052. 

[53] R. Emparan, A. Fabbri, and N. Kaloper, JHEP 08, 043 (2002), hep-th/0206155. 

[54] T. Tanaka, Prog. Theor. Phys. Suppl. 148, 307 (2003), gr-qc/0203082. 

[55] S. B. Giddings, E. Katz, and L. Randall, JHEP 03, 023 (2000), hep-th/0002091. 

[56] S. W. Hawking and D. N. Page, Commun. Math. Phys. 87, 577 (1983). 

[57] C. Misner, K. Thorne, and J. Wheeler, Gravitation (Freeman, 1973). 

[58] B. Kol (2002), hep-th/0206220. 



